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The Load Combination Process 


SOME Abi.) Lyie, Ne felts iis) vo)! sus 


Many physical structures are threatened by combinations 
of loads of varying magnitudes from various sources. For 
instance, bridges, piers, dams, and buildings can experience 
llodids fremyvind; snow, sce, tides, earthquakes, and so 
forth. In many instances the total load or stress exferi- 
enced by a structure varies in time in an apparently random 
fashicn. Certain components of the loads vary rather 
Slowly; others occur more nearly as impulses, such as those 
associated with winds or earthquakes. The problem ~1s to 
design a structure to withstand the superposition of random 
loais- from many sovrces with at least an approximately 
understocd probability. 

The purpose of this thesis is to describe and investi- 
gate certain simple but somewhat realistic probabilistic 
loail models for use in design, and perhaps safety assessment 
Of ISTIUCtUEeS. In this thesis we confine attention tc the 
superposition of just two load types; shock loads, and 
constant loads. For example, winds gusts, flash floods, and 
earthquakes have varying magnitudes and have relatively 
short durations in ccmparison to the times between their 
eccurrences. These will be modelled as instantaneous shock 
loads. Cn the other hand, snow, ice, or water accumulation, 
or even the presence of slowly moving vehicles present loads 
that remain nearly ccnstant in time, occasionally changing 
to new levels. These will be modeled as constant loads that 
change infrequently. Theougheut, this anvestigat vonet caw 
be assumed that the effective stress exerted by several 
types of loads acting Simultaneously can be expressed as a 
Jinear Combination c& “the component loads, the load 


components being treated as stochastic processes. 
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The times between changes in magnitude of the constant 
ioad frocess are independently and identically distrikuted 
with distribution function H. The succesSive magnitudes of 
the constant load are independently and identically distrib- 
weed WIth distributicn function EF. Given the constant load 
process, the shock load process iS a compound Poisson 
process. The successive shock load magnitudes are indepen- 
aently and adenticalty distributed with distribution func- 
tion G. Gre condaticnal rate of arrival of shocks, given 
Enat the Magnatude cf the constant load process is xX, is 

MX) - 

Bet X{t) be the magnitude of constant load process at 
time t; and Y(t) be the magnitude of shock load process at 
erme Ves ACE) —x (L)4Y (£) as the superposition .of the two loads 
at time 62 tard M(t) =SUB Z(s), the maximum load combination 
in Mopper See ibucune at<1. 

Retell, — (in ase Oyen (E)joxas} ae, the: first-passage time 
Zor the load combination to exceed a stress level x. Tx 
will represent the time to faiiure of a structure whose 
EELeCuguaeLs i, sand Wicusubjected to aj stress history 
2. (t Wave s0}i 

Gaver and Jacobs (1981) studied the above inodel for the 
€ase in which the distribution H is exponential and the 
shock oad (mate 4 iS a constant independent of, the shock 
ioad process. Related load combination models that have 
been Studied include Peir and Cornell (1973) { Ref. 1] Wen 
(to) 7) [pret 2) ,Pearece and Wen (1983) [Rek. 3]. 

Asy@protic Tesults which appear in “Gaver and Jacobs 
(i9a4). and Jacobs (1984) will “be described in chapter 2 for 
te GiStribution Of TI, aS xX? ce and Ob the tar: Of the 
CiserioueLon or Tt Lon Linate x. 

in chapter 3, aositulataon, program’ tor the load 


combination model will be described. 


(te) 


| hn re 


Figure 1.1 The Load Combination Process. 
In chapter 4, a Simulation study, of the iwse of, the 


asymptotic results of chapter Zia EO approximate the 


distribution of -T, wallobe idescribed,. 
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Ae IEE DISTRUBUTION CF Ty, FOR LARGE X 


This section summarizes some of the results found in 
Gaver and» dacobs’ (1981) | [Ref., 4] forthe, nodel»iny which the 
Gistribution H is exfonential (a) and the shock load arrival 
rate A is constant independent of the constant load process. 

First, wie tcomsader g the) dastributiuones function votre the 


lidki nun “oadveombinationgtotoccur during [ 0, t)5 0 M(t)« 


H(t) =P {i (t) <x} 
x =P{™(t)<x,Tr >t} +P (M(t) <x,T, <t} ree (2-1) 
where T).1ise the apprival,time,.of the first shock. 
A renewal argument yields; 
P{M(t) <x,T, >t}= E> *[ exp (-ut O-y)) Fedy) 
BM (t)Sx,7, <t}= [tre dul* exp(-uvGix-w) Hyct-v) Ecdy) 
where G(x)=1-G(x). 


Next take Laplace transforms with respect to t of (2-1). 


Hg) = So eft Hatt) dt 


SIH SE Aye Tale (IL aI al haa (2-2) 
where 
Bite Tay Wea ee tere gee Fd yhicese ea === --- 455 ----—— (2-3) 
Hence 
Ect i ol eer Se Lal | ces oe ee ee ao 
HAS = 1 -Mx(¥) si) 


hi.seems.to be difficult Loanvert Hy (E) LOE pan youimt er- 
esting. choice. of the distributiorss F (x).\and” G.(x).- However, 
useful information can still be gleaned from equation (2-4). 
BLESty nant EOm the: definition, Of) Try 

ENG) SSeS eS (25) 

Pinne 8 gS [e822 seri Saal aus | a (2-6) 


Her & 20, 26 equation {2—6) to, find that 


qa 


= Mx (0) 
=i tals = — ee SaaS = 

m (x) B(t ) H,(0) / -X Mx(0) (2 ‘F) 

the next nresult is the JIvmitang property for the f1ese 
time the load combination process exceeds a given level x. 
The limiting distribution of Ty is exponential in the sense 
that 

lin Pim(xy) ,>t}=er sa, one a (2-8) 
XIX. x 
where x, =inf {t;F*G(t)=1} 
The ~ prcot- of “this Tesult for “a more’ = generat medel “thar 
includes’ that ~of chapter ‘2 canbe found in’ vacobs €19GH} 


PReE.- 5 5 


Be TLHE TALL-OP THE? DISTRIBUTION OF “TyebOR FiINatTE & 


In this section we will summarize results for the asymp- 
totic behavior of Pi(I, >t). for large stb -and= faxed" fas top 
xor the model in which the distribution H is exponential and 
the shock load arrival rate is constant. More details and 
results for a more aeneral model-( that’’ancludes that of 
chapter 2 can be found’ in VJacobs (1984) . 

Fron the equations (2%) andy (2-5) 

P (Ty >t) = erti[* ecdy e- Hale Het 


+h de (etaire | eer peasy eee (2=9) 
let 
L (as)= {* crdyy re TAH aOrW)5 ds ------- (2-10) 
A = Ae 
L (= HG F (dd) SeGema a e (At Gx g)1t} ------- (2-11) 
A 
aia Sass (2272) 
The renewal equation (2-9) gives; 
play ot) SGit) R/C Gay Hite) Ae ek Beet ee (2-18) 
where g(t)=[*pcdg) e rtuatx-4jt 0 Bes sei (2-14) 
Following “the” development an “feller” (1971: " page 9376) 


[Ref. 6] it will be assumed’ that there “exists a -*mumber X 


such that 
ge [eX MF Ledyy > eee Pe eee (2-15) 


ttsperocteisa unique and Since the distribution of I is 


defective, is positive. 


Let 
Tce Ise 9 8 eee Sl Cale 
gi fe GMOS Ge a ee (2-17) 
Pe(ty>t)= e*PTHY (t) we eee een n--------- (2o18) 
Then 
P¥ (Ty >t)=g#(t)+ [PP (T,>(t-s)}L* (ds) 9 --------- (2-19) 


It is assumed that the assumption for the key renewal 


theorem holds. Aprflication of the key renewal theorem 


yields 
. Kogt ey LE Cee Tae | ye a ee fae ae eS 
Lime Baty >t) = ra 9 (t)at (2-20) 


G4 (tyat= fr (ay) [Met oMaowIE at 


ifs ° 
/ 


_ (ee Memento) we NO ar Ce Se : 
=! F(dy) ATM &le-$) - K(X) ng 


Rewriting the equation (2-15), the equation determining 


K is 
Be ee | ed eee SEAS dlc 
= : 
x n»~ 
= cd eS CU, = - 
[" ecdy oo eeerae (2-22) 
and 
ee pele ne Ne OE ode 
= (* J ometh 2, Ue Sheena See ee SSeS 2253 
J pay [A+ a@er-g) - Kx) ]* 
The key renewal theorem implies; 
Lim e*%#t>(T, >t) = — : Tag a Raat -- (2-24) 
sr No FCdy) ——_>__________ 
Exanples: [ A+ UM G@lX-4) - KU) 


ie eappeats. tolpbey dititicult: feoysolve; for kK analytically 
in general. We will discuss the numerical computation: of K 
fotts dm wo) | examples: In both examples, A=A=1 and the 
dtatrabutvon Or the shock and constant load magnitudes Firand 
G @isel ofthe form: 

EF (i= eo" 


GGj= e-o% 


In example 1, a=b=1 and in example 2, a/b=2. 


To. compute the & for F(x] cs ee (ae - 
et 
A 
Lig: eee SN Se ed es” . 
EGS EA) eemege ae (2-28) 


where A=A=1. 
Her model” 1, 
f (x)= TRE {U2 07)" [Uen 1 404) = Lens KOH e~ py ==, (2-6 
for K<lt+e™* andi 305 
FroriK=T;rt @)\ has’ a Gingullarspoint;: so hy usangetne 
lL HOPITAM REE fRelee 7 | 
Zam f (x)= Lim g(x)/h'(x) 
where g() is the numerator and h(x) is the denominator 


Of the Funectson = (7: 


For K=1, 
fj) —— Ke —=} ---- (2-27) 
For model 2, 
-bx 2 obey -bx 
£ (kK) = 74 es [bias x+e™%)* (= a Ne - ie —x)UIH e)) 
£00 (ace eo) ae ie pe Sylp Aha LY Se See (2-28) 


for x <1+e™ and K #1. 
For X=1 £() has also singular point, so by using tre same 


method as used in model 1; 


HORIG—=a 

£(1)= aS (e- ™) == (2-29) 
Differentiation of the equation (2-25) gives 

fig afte aa ae 


[ A+ am eb&s -J> 
Fouation (2-30) ws igneater than 0) kor O<x<1+e™% and x1. So 
the function £ (x). is as monotone. increasing?’ funetaon am the 
interval! 0<*<1+e% 5. Using, thisrmesult!) (¥ we! computelth ork 
from the equation (2-22) by applying the bisection method 
[Ref. 8]. 
Using the computed kK we compute the constant y¥* such that; 


# = i cd Se ee 
y Je FCd#) CL A+ MG -4) -KJZ 


is the constant of equatwven (2-24). 
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Fer case a=b=1, 
: oe = -x)_2e% = Rhea 
OF = Gage {[(i-xte%) -2e a -K te exe) 


-[e-* (2-x)-2€7* In (OE 2-) — i ae S| -- (2-31) 
For case a/b=2, 


-4¢400 ~x 
———_ titer Paet)s Pee (er (ee) 


# 
ghar. 


=ib I -3b 
$3I€ 24 “(den (ie K+ 07>*) Aon ((I-K Pee enpehe= at U-KeMe™ 1- ries ac #2) 


III. COMPUTER SIMULATION 
This section describes the computer Simulation model. 

We denote thejeconstantiloadl qusiveat erates haan), the. shock 

load arrival‘rate by ,~,, “the dastribution cf constant Load 

Maghitude by F(x), and the distribution of shock load magni- 

tude by G(x). The simulation model will be used to study 

the quality of the approxamations of “the “distribution of 7 

by the asymptotic results described in chapter 2. 

The following describes the simulation algoritkhn. The 
stress ievels are X, <x,<X,~.---.-.--<X,- 

estate =O ie =(0.4 

2, Set V=0, generate a constant load magnitude,C, and 
duration pnterval Te. Find the largest: andex “1>N, 
SuUChS chae ex SC 
Put N, equal sto that andex and 
set 

ty =f Gor, <uShy, 

Set No oNy ie 

3, Generate a shock load magnitude,S, and a time until 
the shock occurs, Ts. 

4, If the shock interval does not exceed the excess life 
of the ‘constant load duration, Tc, find the, Langest 
index’ 1>N, "such ) thatex,.<ctS. 

Let N, be that index. 
Set TyeHPTtVtls for Nealon. 
set No =N,.,V-Vtas,, [c-Te-i1s. (Go Lo sstep -3. 

5, if the shock interval exceeds the excess life time 

Of the constant sload durataton, 


set E=T+tl¢, go toustep Z- 
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One realization of the simulation is completed when one 
Ty 1S €COnmputed for cach Level x }1Sisn. ‘Back Simulation 


consisted of 5000 replications. Further details are 


ee oe ee ee 


| 


Fagure 3.1 Generating T. 


included in chapter 5. 
The key subprograms are following. Other programs will 


be explained by reading from the output results. 


RANDO 3; generate the Tx. 

CCMPUT ; compute the P(Tx=0), true P(Tx=0), sample nean, 
true mean, standard deviation, coefficient of 
variance. 

Sample probability 
Ede -Nalwesparapae var is obtained from 
the sample. 

2rue ypro habad ety, 
E (Tx =0) =1-P (Tx>0) =1-P (M(0) $x) =1- frF(€y) =F (t) is 
obtained from the equation (2-5). 

The sample mean 


ie A 
Xx =e, exe: an 


uly 


the true mean E(T,) is obtained in special 
cases from the equation (2-7). 
See details in Gaver and Jacobs (1981). 

Sample standard deviation 


~ nN 
(\ = er lle is obtained from the sample, 
i= I : 


Sample coefficient of variance 
COW) c= a is obtained from the sample 
mean and sample standard deviation. 

CCMP ; compute tke sample guantiles, and the guantiles 
cof an exponential distribution with mean E(Tyx). 

A sample quantile is 

=(dxSanple size) th order statistic of sample and 
bt = ~-ECTy) Im( 1-4) is the guantile of exponential 
distribution with mean E(Ty,). 

CCM ; compute the quantile of the sample conditional 
distribution given Tx>0 and the exponential distrib- 
ution with mean E(T,)/P(Tx>0). 

The sample quantile of the conditional distribut- 
ion 
9 e=[ 4x (number OL positive Ty) Ji hvordersstarrstac 
is obtained for sample cf positive Ty. 

The conditional quantile of the exponential dist- 
ribution with mean E(Tx)/P(Tx>0) is 


| EGR 
bac = PC 7% >0) Ln (I5dc) - 


SK 3°compute the *K such en 
K={ K(x): v=]? cedy) 
bisection method. 


BEM ae Y)- we by applying the 
INT ; compute the constant y* such that ; 


Y?=oCtdy =e 
‘ J #) CA+UGU-¥) -K 2 
by substatuting the KX Hound=in Sk. 
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The following cases were simulated. 


arser i burron case EN AAC) Fax) G (x) 
a,éxponential hie 1 1 1-e” 1-6” 
(a=b=1) Dy 1 2 " " 
3, 2 1 WwW W 
b,é€xponential le 1 1 T=er4” 1-e%% 
(a=2, b=1) Dy 1 Y W W 
ee BY 1 n W 
c,expenential a 1 VAG 1-e™ 1-e™% 
with varying Dy 1 ea Hu Wy 
arrival rate 3, 1 ec M " 
d,Fareto case di 1 2 1-=y 1-e%* 
2¢ 1 2 lagen tae 
/ 
<1 1 : Uae I~ TR 


These different settings will be referred to as case (A-1), 


eo eee Daly B=) en iKB=o), (Hl), (C2), (C3), 
P-ljyte-2), (D=s) Trespectively throughout this paper. For 
€ach above case, we ccmpare the asymptotic result (2-8) and 


the experimental results. 
The next setting is to examine the exponential approxi- 


Bae von Or etre Stari cf the distribution of T, for finite x. 
the cases of identical exponential distributions, 
F(x) =G (x)= e* and different exponential distributions, 


F (x)=e™ G(x)=e* a/b=2 


where 
were considered. 


the E£ollowing conditions are considered. 


Fh 


=e 
Is = = 


Bix) —C¢ (x) —e.- 


F (x) =e"** 
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,G (x) =e™* 


These different settings will be referred to as case (£) and 
case (F) respectively. For each case, compute the #&, and 
constant 7 and compare the guantiles of the asymptotic 


results (2-24) to the simulated data for each level of x. 
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IV. ANALYSIS OF THE SIMULATION RESULTS 

NS 1 nda Cared A ant chapter, 2) tthe .xact? distributivon’ jot) Ty 
ised EFI Cults to) obtain in * general.” |) ini this chapter si mula= 
ton, wild | be-usied) tor study sthe accuracy of » the).approxima- 
Eaions "0% sbhe dastribution Pot 1 Ty i wby the exponential 
distributions suggested by the asymptotic results of chapter 


ie 


Ae ebHe  DESTRIBUIEON: (CFO Tx) FOR, INCREASENG .X 


@gne Itmiating result (2-8) inducateduthat the tistEs bu— 
tion of Tx approaches the exponential as level x increases. 
tits) SLeoul et SUGgGeSts, that the distribution of Ty can be 
appEeoximateds Dylanjyexponential dvstribution Jat! least for 
large x. In this section this exponential approximation will 


ke studied via a simulation. 


— = = SS Se SS a i a iO Sh SS ee er 


ae = SS SH ee SS SS ES SS a 


i Trice wrest COLLECtIonN, Of “three, models: to be 
considered F (x)= (x) =e. The times between constant load 
changes are exponential with parameterAand the shock load 
aEbivdal Eater “25 a2 constant, independentes o£) they= constant 
load process. 

For these cases it iS possible to derive an analyt- 


ical-.cxpression for &(Tx'). 


E Uf y= se ey -UXe% +Me™ In ==) | ee ee ee ee (4-1) 
ines aa ee 3 
~ n» A+uM e-* 


ie Ehews wa beLes, 4a P {Ty » =0) refers:,to,the simulated 


Samp Veasprebabida ty, pthadtyutT ~=0s> ) X= BAR, | Gives; they ,»samudated 


sample mean of 'Ty ¢ “ST=DEV “US the: simulated? (standard 
deviation-okrr Ty 7.” SVAR (XB) is the variance of the sample 
mean; COEF-V is the sample standard deviation divided by the 
sample mean. T P(Txy=0) is the theoreticai probability that 
Tx=0, F(x); T x-bar is the theoretical mean of Tx computed 
from the equation (4-1). 

The simulated coefficient of variation in all three 
cases decreases as the level x increases becoming quite 
close to the theoretical exponential value of 1 when x=5. 

To compare the simulated distribution of Tx to the 
approximating exponential distribution, quantiles from the 
Simulated data were computed and compared to the corre- 
sponcing quantiles cf an exponential distribution having 
theoretical mean (4-1). The Simulated quantiles were 
computed as described in chapter 3. The exponential 
.d&-quantile is given Ly 

g, =~a-ln (1-4) 
where a is the mean of the exponential. 

For each level of x, the Pirst. row Gn the tables 
(A. 18,2 (A. 2) GL (AC BYE fOr thes Gquantales oft heidasteabutien of 
Tx gives the quantiles of the simulated data and the second 
row gives the corresponding exponential quantiles for an 
exponential distribution having theoretical mean (4-1). 

Cne way the distribution of 2x dvtfiers from an) iexpo— 


nential is that 4 has an-atom at 0. in parevcuvan, 


PAT =O} =F (x) —ean 
for the case considered here. 

Quantiles were used to compare the simulated condi- 
tional distribution of Ty givenefy->0 “to ant “expomwentaa 
Gistribution haying mean B'(Da)/P (2a O)a . “Phe simulated gquan- 
tiles for the conditional “distrrbution |were Computed sais 
Gescribed in‘“chapter 3.- These* quantrles* were’ compared te 
those of an exponential distribution having theoretical nean 
B/E Cire Or. 


Ze 


in the pt ables Aad). Amis Ze cand oy(Mar3)- o eo sthe »quan- 
tikesyot the .condiitvonal sdistributvon,£o0n each level »x;-~- the 
first row shows the simulated quantile and the second row 
the corresponding approximating exponential guantile. 

The quantiles of the simulated conditional distribnu- 
tion are much closer to their exponential approximation than 
tie Guantztes) fer the cunconditional distributzon » sto«their 
exponential approximation. However the quantiles for the 
mcondi tional distribution get closer™™“to “those of their 
approximating exponential as x increases. 

Gonpasa ma gtahles o(Mal),)), (Ao2)i, -ands(A Sy at appears 
Boat. inereasing.gthe sarrivals cate of jshockh) jor the <rate of 
change of the constant load magnitude decreases the quan- 
tiles. The change of arrival rate of shocks appears to have 


the «greater veffect. 


In the next three cases, Ehetdarsitrebutawvon #yor wthe 
shock ioad magnitudes is exponential with parameter 1 and 
the distribution of the constant load magnitudes is exponen- 


tial with mean 0.5. The other model assumptions in tables 


(Bil hy (Ba2)i, (BJS jew commesp ond) sto» sthosesimitables p(As ln; 
(ian2) grand 4 (h.:3))pe 
As before, for each level x, the phest «row Ain «the 


guantile table gives the simulated guantile and the second 
row gives the approximating exponential quantile using the 
exponential distribution with theoretical mean. Comparing 
the quantiles .o£ Tisingtable», (Bek) swarth)(A. 1) (respectively 
(242 yn daa Ss (a2 jon and» (B.3) Wisttion (eS )y)g) sakes seen that 
decreasing the mean constant load magnitudes has increased 
tie iqudntiles. 4 rusthensz) it Jappears thatthe convergence of 
the distri bution a0 fel) to gam 4 fexponenitial vis daster, in the 


case of smaller mean constant load magnitude. 


a3 


The exponential approximation to the conditional 
distribution of Ty given T,>0 has theoretical mean 
E(2Ly PAELLR D0) s Once again for small levels x,it appears 
that the exponential approximation to tthe “Conds trtonat 
distribution of is better" “than that forthe’ unconditional 


one. 


load Magnitude 


In the next 3 cases, constant load magnitudes and 
shock load magnitudes are exponential with mean 1. Constant 
ioad magnitudes change according to a Poisson process with 
rate 1. “Given the ‘constant lload ‘magnitudevat time’ t*isee; 
the prokability a shock loadewHld dareives tne teeters 
interval of t,,bth 9 ses. Ac) hn +o) Ch). 


In case (C=1) the conditional shock load arrival 


rate is Al (cy =C< Comparison with«table ({(A.1) tndicates that 
conditional shock arrival rete yic) “=€ has the sefftectios 
decreasing the quantiles of fT,. Further the exponential 


approximation “to ‘the quantiles tof Ty Vis "net tas igood- ras tras 
case’ "A—1)\. 
The approximating exponential distribution to the 
conditional distribution of Tx given T,>0 has a mean of 
E(Ee Ze (1x >0).< The quantiles of the exponential approxima- 
tion to the guantiles of the conditional distribution appear 
to be closer than those for the unconditional distrikution. 
In case ©(C=2)-- shock loads “arrive” jwith condi: tronayz 
appival@rate satcy! l=expi(ce Comparing the table (C.1) and 
(C.2), it is seen that the ‘quantites -of the Case w(cpe—e are 
less than those for the case u(C)=exp(C). This is because 
interarrival times tends to be larger for the case 
MC PP =expicy: Comparing ‘tabies 7 (An) wand "(Coz) indicates 
that the exponential approximation to the quantiles of ‘x is 


not as good as in the case M=1. 
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ine Casey (Eze) ““the sconditaiondlyishockvarrival ‘rate irs 

INES Sexmpi(=C)< Comparang, tablesed(Cs2)y “and? (C23) yo indicates 
that the guantiles for the case 

A(C)' =exp(-C) are less than those for the case y(C)' =exp (C) 

and furthermore the exponential approximation of Tx appears 

tePle+elose fer “the "case ?(C=3) than that of "cor the case 

(e= 2). In all of the cases the simulated mean was used as 


the parameter in the exponential approximations. 


In the next three cases, the times between constant 
load changes of magnitude are exponential with mean 1. Shock 
loads arrive according 40%a POLSSON “process with “rate,2. “In 
case Ma) ) thes sconctaut’ load smagnitudesshave «a «Pareto 
distribution with parameter <=1 and the shock load magni- 
tudes are exponentiaily distributed with mean 1. Comparison 
Weeh Seabie.~ (A<«2) andicates, that, the Pareto, constant; load 
magnitude has the effect of decreasing the quantiles of 1f,. 
Further the exponential approximation to the quantiles of ] 
is not good in the Pareto case; (the approximating exponen- 
tial has a mean of tke simulated sample mean). 

the approximating exponential distribution of %& 
given Tt >0 to the conditional distribution has a mean of 
EAT sai S Pi igs SD) 2 Once again the guantiles of exponential 
approximation to the guantiles of the conditional distribu- 
tion appear to be closer than those for the unconditional 
distribution. 

Ee Case s4(D=2) the constant load magnitudes have a 
Pareto distribution with parameter A=2. Comparing takles 
(Pe tjee duds (>. 2) eatsis seen that,-thes quantiles of, the sinu- 
lated distribution cf Tx, for the case &A€=2 are larger than 
those for the case cdA=1. Comparing tables (D.2) and (A.2) 
indicate that the exponential approximation to the quantiles 
of Ty 25 MOte) dS, good; ;£OE4. Lhe.Pareto,. case ase tors the 


exponential case. 


DD 


In’ ‘case*((D=s) the constant load magnitudes have 
Pareto distribution with parameter «=1 and the shock load 
Magnitudes have Pareto distribution with d=1. In all three 
cases the simulated mean was used in the approximating quan- 


ties 


B. -THE “EXPONENTIAL TALL OF THE “DISTRIBUTION OF ‘Ty,’ FOR 
FINDTTE <x 


In this section sSimulaticn will be used to study the 
exponential approximation 


suggested by the asymptotic result 
Tin Cte) = 4+ 


This’ iS an approxinatacone for P\(ly >t) MOE £i xed Elnate) xe 
Should become more accurate as t becomes large. 

Two cases are considered. In both cases shock loads 
arrive © according = to’ a’ Poisson process: with) rate 12) end 
constant load magnitudes change at the times of arrival of a 
Poisson process with rate 1. The shock load magnitudes have 
exponential’ distribution’ with medn~ 129% in® ? case? (a), the 
constant loai magnitude distribution iS exponential with 
mean 13; In case (F) it is exponential with mean 0.5. 

Descraption~ o£ ‘the’ computation of XK and t# for these two 
cases can be found in chapters 2,3. It can be seen from the 
computed vaiue of x(x) ana ¥ (x) appearing in tables, (E24) 
and® (Fol) ¢ that 'K(x) Mane approaching Phe valve 7b (in (where 
E(Ty) is the theoretical mean) as x becomes larger and for 
the cases computed K(X) SF -o OFureE hier | ¥(x) is approaching 
Iwas’ xk becomes” largex. 

-To study the exponential approximation to the: distribu- 
tion of TI , Gguantiles from the Simulated data were computed. 
These can be-found in tables “E22) and’ (F22). Ferredch Veved 
x, the first row presents the simulated quantile; the second 
row gives the approximating d&-gquantile computed by 


q,=eckt-in{P(i-d)}] 000 cee ---------------------- (Toles 


2d KON) 


The third row shows the approximating %-quantile computed by 
See eae i a eae ane == = = -=—- (eel 2) 

where the theoretical mean is used for E(T,). 

€Ceonmpariangeehe gwantiles, ian tables (E.2) and (E.2).it-can 
ke seen that the approximating quantiles (LZ) are close to 
the simulated ones even for <= .40. The approximating guan- 
tiles computed in (qf) do less well but improve as x becomes 
larger. The two approximating guantiles approach one ancther 
as x becomes larger, as expected. Approximating quantile (Zz) 
has the advantage, however, of being eaSier to compute. 

hapless (G. ij. presents 90-7", “comfidence intervals for 
da-qguantiles with 420.9 for case (E). These confidence inter- 
vals were computed using a large sample approximation, see 
Conover “(1930s page 1171-1.14)-of Re£.-9 J. 
Comparison with table (E.2) suggests that the approximating 
guantiles approximates the simulated ones well for these 


values. 


Cc. . CONCLUSIONS 


Approximating guantile I is always better than approxi- 
mating guantiley. Approximation L approximates q, well for 
as Small as 0.49. However J is more difficult to compute 
that E+. 

The “pereoEmance OL approximating ~ quantile 2 | can. ‘be 
improved by changing it to 


= ECTx ) 


eR aaa 


"4 PC Tx >o) 
and using it to approximate the quantiles of the conditional 
dastEi bution, o£. ty guven, 1%>0;. for. the models considered 


P(™ =0) is easy to ccmpute being egual to F(x). 
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Lc=constant load magnitude 
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Ts=shock load interval 

Ipc=ccunting the number of repetition 


L=counting the maximum level of X; 
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